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Abstract —Target tracking faces the challenge in coping with 
large volumes of data which requires efficient methods for real 
time applications. The complexity considered in this paper is 
when there is a large number of measurements which are 
required to be processed at each time step. Sequential Markov 
chain Monte Carlo (MCMC) has been shown to be a promising 
approach to target tracking in complex environments, especially 
when dealing with clutter. However, a large number of mea¬ 
surements usually results in large processing requirements. This 
paper goes beyond the current state-of-the-art and presents a 
novel Sequential MCMC approach that can overcome this chal¬ 
lenge through adaptively subsampling the set of measurements. 
Instead of using the whole large volume of available data, the 
proposed algorithm performs a trade off between the number 
of measurements to be used and the desired accuracy of the 
estimates to be obtained in the presence of clutter. We show 
results with large improvements in processing time, more than 
40 % with a negligible loss in tracking performance, compared 
with the solution without subsampling. 

I. Introduction 

Flooded with data, richly provided by electronic sensors, 
the current monitoring systems face the problem of being able 
to process the data and monitor the phenomenon of interest 
at the same time. In this paper we consider the problem of 
target tracking in large volumes of data. There is a wealth of 
algorithms that can provide sequential estimation of the states 
of the target, e.g. for details see HI. 12- In a Bayesian frame¬ 
work, the posterior distribution can be iteratively computed. 
However, analytically this can be achieved only when the state 
space model is linear and perturbed by a Gaussian noise. In this 
case the solution is referred to as the Kalman Filter. There are 
a large number of techniques which overcome the limitations 
of the Kalman filter based on the sequential Monte Carlo 
(SMC) methodology. The seminal work on SMC in target 
tracking was presented in go which was also referred to as 
the bootstrap particle filter (PF). The bootstrap PF and many 
variants thereof, broadly referred to as PFs, are commonly 
favoured techniques in a wide variety of applications due 
to the filters ability to handle non-linear state space models 
and/or state space models perturbed by non-Gaussian noise. 
However, the PF is not void of challenges. Some of the 
difficulties faced by PFs includes weight degeneracy and 


sample impoverishment. Although there are variants of the PF 
which have been proposed to alleviate these issues ID, 0, the 
PF is still susceptible to degeneracy, and these difficulties are 
more profound when tracking complex systems. 

Markov chain Monte Carlo (MCMC) techniques are a 
powerful set of algorithms for sampling from a probabil¬ 
ity distribution. MCMC tecnhiques, such as the Metropolis 
Hastings (MH) algorithm, have been predominantly used in 
applications requiring static inference 0. Recently there has 
been considerable interest in extending these techniques to 
sequentially updating the posterior distribution (7), (8). Se¬ 
quential MCMC has shown promising results for complex 
systems. The largest hindrance being long processing times 
which could limit usage in applications required to run in 
real time. There have also been several algorithms ED, EDI, 
02 which have been proposed to help reduce computational 
complexities when performing static inference with MCMC 
techniques on large datasets. 

In this paper we propose a novel technique which results in 
an efficient sequential MCMC algorithm when applied in com¬ 
plex systems consisting of a large number of measurements. 
This is achieved through the combination of sequential infer¬ 
ence and adaptive subsampling of the measurements at each 
time step. We show how the proposed adaptive subsampling 
sequential MCMC algorithm can be applied to target tracking 
and illustrate the computational savings it affords. 

II. Problem Formulation 

Target tracking of a complex system can be considered as 
sequential state estimation with multiple measurements. This 
can be achieved in a Bayesian framework by sequentially 
computing the filtering posterior distribution p(xk\zi-.k ) where 
Xk € 11™“ is the state vector at time tk with k = 1,..., T £ N, 
and z\-k = {Z |,..., z, represents all the measurements 
received up till time tk- The measurements received at each 
time tk are represented by a set Zk = {z \,..., z^f k }, where 
Mk is the total number of measurements and z\ £ R n *. The 
filtering posterior distribution can be recursively updated based 



on 


p{xk\z\-.k) ck J p(zf i \xk')p(xt c \xj i —i)p(xk—i\zi : k—i^dxt £ —i, 

(1) 

where p(zk\xk) is referred to as the likelihood probability 
density function (pdf), and p(xk\xk~i ) is referred to as the 
state transition pdf. An analytical solution to (Hi is typically 
intractable when the state space model is characterised by non- 
linearities and/or non-Gaussian noise. 

A. Sequential Markov Chain Monte Carlo 

MCMC methods work by constructing a Markov chain 
with a desired distribution as the equilibrium distribution. A 
common MCMC technique used to obtain samples from the 
equilibrium distribution, n(x), is the MH algorithm. This is 
achieved by first generating a sample from a known proposal 
distribution x* ~ g(-|a: m ~ 1 ). The proposed sample is ac¬ 
cepted as the current state of the chain, x m , if the following 
condition is satisfied 

ttQ v*)q{x™-'\x*) 
Tr(x m ~ 1 )q(x*\x m ~ 1 ) ’ 

where u represents a sample from a uniform random variable 
u ~ C/[ 0 , 1 ] - Using Bayes’ rule and assuming that there are 
M conditionally independent measurements, z 1 , results in the 
further expansion of this expression 

p{x*)q(x m ~ X \x*) p{z*\x*) 

p{x m ~ 1 )q{x*\x m ~ 1 ) p(z l \x m ~ 1 ) 

The previous state of the chain is stored as the current state, 
x m = x m ~ 1 , when the proposed sample does not meet this 
criterion. We further manipulate this expression into a form 
with the likelihood isolated: 


as the equilibrium distribution in order to help alleviate 
the high computational demand. In a similar fashion, an 
approximation for the joint filtering posterior distribution 
can be obtained through MCMC methods by representing 
p(xk-i\zi-k~i) with a set of unweighted particles. This ap¬ 
proach has the advantage of avoiding the direct Monte Carlo 
computation of the predictive posterior density. Furthermore, 
the approximation can be marginalised to obtain the filtering 
posterior distribution of interest. 

More specifically, at each time step, the particles are updated 
with a MH joint draw for Xk and Xk-i, followed by an 
individual MH draw for Xk . The second step, referred to as the 
refinement step, is introduced to aid in the mixing of the chain. 
An appropriate burn in period, N}, urn , was also introduced to 
minimize the effect of the initial values of the Markov chain. 
This results in the definition of the total number of MCMC 
iterations at each time step, N = N p + Nt, urn - This approach 
is highlighted by Algorithm Q] and is referred to as standard 
sequential MCMC. This approach showed promising results 
in a multi-target environment but is still susceptible to high 
computational complexity when a substantially large amount 
of measurements are required to be processed. 

B. Adaptive Subsampling 

In standard sequential MCMC, it is required to perform 
2 NMk calculations of the likelihood at each time step. This 
is highlighted in the computation of the log likelihood ratio, 
Af fc (-) and Af fe (-), in Algorithm [T| When Mk is very large, 
the log likelihood ratio becomes the most computationally 
expensive step of the algorithm. To reduce the computational 
complexity, we introduce a Monte Carlo (MC) approximation 
for the log likelihood ratio: 


log 


p(x m 1 )q(x*\x m *) 
p{x*)q(x m ~ 1 \x*) 


M 

< E log 

2 = 1 


p(z l \x*) 
p(z l \x rn ~ 1 ) 


(4) 


In 0| it was proposed to use MCMC methods, specifically 
the MH algorithm, to target the filtering posterior distribution 
in 0 as the equilibrium distribution. This allows for the 
iterative update of an approximation of the filtering posterior 
distribution by representing p(xk-i\zi : k-i) with a set of 
unweighted particles. 


iV„ 


p(x k -i\z 1:k -i) 


N, 


■E 5 ^- i-4-i 

3 = 1 


(5) 


where N p is the number of particles and (j) the particle index. 
This technique was shown to work well in state space models 
containing a high number of dimensions when compared 
to techniques relying on importance sampling, however, this 
direct approach may result in a high computational expense. 

It was proposed in (8) to consider targeting the joint filtering 
posterior distribution of Xk and x^-i 


p(x k , Xk-x\Zl,k) oc p(Zk\Xk)p(Xk\Xk-l)p{Xk-l\zi:k-l), 

( 6 ) 
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(7) 


where the set z* k = {z ^’*,..., z k m ' k '*} is drawn uniformly 
without replacement from the original set of Mk measure¬ 
ments. 

The difficulty which arises is in selecting a minimum value 
for S mt k that results in a set of subsampled measurements 
that contain enough information to make the correct decision 
in the MH step. To overcome this difficulty in standard 
MCMC for static inference, the authors in [ 10] proposed to use 
concentration inequalities which provide a probabilistic bound 
on how functions of independent random variables deviate 
from their expectation. In this case, the independent random 
variables are the log likelihood ratio terms. Thus, it is possible 
to obtain a bound on the deviation of the MC approximation 
in (7} from the complete log likelihood ratio: 

P(|Af^ (x™-\x* k )- Af* (x™-\ x%)\< c Sm , k )> l-5 Sm , k 

( 8 ) 

where Ss m k > 0, and cg. m k is dependent on which inequality 
is used. There are several inequalities which could be used, in 















Algorithm 1 Sequential Markov Chain Monte Carlo 
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Initialize particle set: {aio^ljit 

for k = do 

for m = do 

Joint Draw 

Propose {xl.xl^} ~ q x (x fc ,a; fc _i|a 
Compute ipi(u : x * k , x k _ x , x™~ % , x™Zi ) 




n—1 \ 
z-1 ) 


= I7T lo g 


p(**l x iE-i)p( a! fc-il*i=*-i) 


7: Compute A.^ k {x* k ,x™ x ) 


= — V M ‘ log 
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ifAf^K^r 1 ) 

> 1 pl(u, x 

{ X k l ^ X 'k-l} = { X ki X k -ll 

else 

end if 

Refinement 

Propose {at^} ~ q 2 (x k \x^,x^_ x ) 

Compute ip2(u, x* k ,x r k ,x r k _ x ) 

P( X k L \ X k‘-l)<l2(xt\xF ,=c]"-l) 
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log 
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16: Compute Af fe (xjj\ a£) = 4- E^i log 


p( z !U a 


.PC**!®™), 


if AjA (x k ,x™) > tp 2 (u,x* k ,x^,x^_ x ) then 

™m _ _,* 
x k ~ x k 

end if 

if to > Nb urr , then 

~ t { m ~ Nburn) _ „m 

x k ~~ x k 


end if 
end for 
end for 


P{ x k\Zl:k) = 4~ E^l S ( X k - X { k :) ) 
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this paper we make use of the empirical Bernstein inequality 
tm ed, which results in: 


CS m , = 


2Vs m , k log(3/6 Smik ) , 3R k \og(3/6 Sm , k ) 


(9) 


where Vs m k represents the sample variance of the log likeli¬ 
hood ratio, and R k is the range given by 


Rk = max 

1 <i<M k 


log 
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Looking back at the standard sequential MCMC approach, 
we find that the joint draw is accepted based on the con¬ 
dition A^(a£,®£**) > Mu,x* k ,x* k _ 1 ,x™~ 1 1 x™jr i 1 )- It 
is required to relate this expression in terms of the MC 
approximation of 0. Since the MC approximation is bounded. 


we can state that it is not possible to make a decision 
when the value of ^i(u,x k , x k _ x ,x™~ x , x™~ x ) falls within 
the region specified by the bound. Thus it is required that 
|Af m ’ fc (x^-\x* k ) - Vh(u, x* k , xi_ x ,x?-\x^)\ > c Sm<k 
in order to be able to make a decision, with probability at least 


This forms the underlying principal for the creation of a 
stopping rule na, on. Let 5 G (0,1) be a user specified 
input parameter. The idea is to sequentially increase the size of 
S m k while at the same time checking if the stopping criterion, 

\4 m '\*'r\ x i)-'k(u,xt. 




u k -1 


> C S„ 


is met. If the stopping criterion is never met, then this will 
result in S' TO ,fc = Mk, he requiring the evaluation of all 
the measurements. Selecting Ss m k = 4 1 » 6 results in 

Es m „>i 6 s m , k < S. The event 


pSr, 


£ = 


fl < c Sm , k ] 

^ m,k ^ 1 
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thus holds with probability at least 1 — S by a union bound 
argument. 

This iterative procedure allows for an adaptive size of the 
number of measurements required to be evaluated. However, 
there is cause for concern with the definition of the stopping 
rule. That is the fact that the range, R k , used in the calcu¬ 
lation of 0, is dependent on the log likelihood for all M k 
measurements. Calculating this range would thus inherently 
require at least the same number of calculations as in the 
standard sequential MCMC approach. In certain applications it 
may be possible to obtain an expression for the range which is 
independent of the measurements, however, this is not the case 
for the current application of interest. In order to overcome the 
computational complexity of the calculation of the range, and 
to reduce the sample variance Vg m k in the bound, a control 
variate has been introduced in referred to as a proxy: 


Pi(x k l ,4)~ lo g 
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Thus the MC approximation in 0 is augmented into 
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It is required to amend the MH acceptance accordingly to take 
the inclusion of the proxy into account. 

We propose using a first order Taylor series as an approxi¬ 
mation for the log likelihood, £,,(x) = logp(z*|a:), given as 

h{x) = £i{x + ) + (V£i)E ■ (x - x + ), (14) 

where (Vf*) a .+ represents the gradient of Ifix) evaluated at 
x + . This results in the following form of the proxy 

Pi{ x ™^\ x t) = h{ x t)-h{ x ™~ 1 ), 

= (V^ + -(®iE-*r 1 )- (15) 








































With the inclusion of the proxy, the range, R k , is now 
computed as, 


R k = max 


P( 414) 


fcl x k) 

l<i<M k [p(z i k \x^- 1 ) 


l<i<M k Yp(z\\x™ ) 
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, x t) 
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We can derive an upper bound for the range, RJ;, i.e where 
k > Rk, which can be computed efficiently 

444) 


R? = 2 max 
l<z<M fc 


log 


pKK" 1 ) 


— pi\X k , X k ) 


2 max {k(4)-^(4T ^l} 


i<M k 

=2 max {|B fc (x™ _1 ) - B fc (Xfc)|} 


(17) 


where B k {x) = ^(a;) — £j(a;) is the remainder of the Taylor 
approximation. The Taylor-Lagrange inequality gives us an 
upper bound on the remainder term. More specifically, if 
|j(n+ 1 )( a ,)| < y, then \B k (x)\ < Y ^~+i)! —, where in our 
case n + 1 = 2. Upper bounding the Taylor remainder finally 
results in the following upper bound on the range 

R b = 2\\B k (x™- 1 )\ + \B k (xl)\\, (18) 

which is dependent on the maximum of the Hessian of the log 
likelihood, Y. The complete adaptive subsampling sequential 
MCMC approach is illustrated by Algorithms 0 and [3] 


III. Application to Target Tracking in Complex 
Systems 


A. Target and Sensor Modelling 

In this application the state vector consists of the posi¬ 
tion and velocity of the target in a two dimensional space, 
x k = [x k ,y k ,x k ,y k ] T ■ The target motion prediction is 
performed according to the near constant velocity model. This 
results in the state transition density having the form 


p{x k t) —Af {x k | A k x k —\, Qjk ), 


(19) 


where Afi(-) represents the normal distribution, and matri- 

' I 2 T s I 2 ' 


ces Ak and Q k are defined as A k = 


0 2 


and 


Qk = 4 


, where T s = t k -t k - 1 . 


(T s 3 /3)J 2 (T s 2 /2)J 2 

_ (T?/2)h T s I 2 

In this application, the total number of measurements re¬ 
ceived is given by M k = M k + M k , where M k represents 
the number of target measurements, and M k represents the 
number of clutter measurements. The number of target and 
clutter measurements are Poisson distributed with mean Ax 
and Ac respectively. The likelihood density thus takes the 
form IT5i : 


M fe 

p(z k \x k ) (X IJ X x px(z k \x k ) + A cpc{z l k ), (20) 

2=1 


Algorithm 2 Adaptive Subsampling Sequential Markov Chain 
Monte Carlo 

1: Initialize particle set: 

2: Determine initial proxy parameters. (See Section UlI-BI for 
more details.) 

3: for k = 1 ,...,T do 
4: for m = 1 do 

5: Update proxy parameters. (See Section lHI-BI for more 

details.) 

Joint Draw 
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Propose {tc*,®*^} ~ q 1 (x k ,x k -i\^ k 
Compute ipi(u, x k , x k _-y, x™ -1 , x™- 1 


™"1-1 

Xu ) 


h k -1 


= A l °z 


p(K\ x k-i)p( x k-i\ z r.k-i) 


n 1 ) 


fe_1 

Compute A^ix^x™- 1 ) and {pi(x^~ l ,x* k )}^L\ 

with the routine described by Algorithm [2 

>Mu,xl,x* k _ v x^ \x%_l) 

-itt sf=i Pit** 1-1 .*:) then 

{ x1 k i x 7 k-\} = i X k’ X k -ll 

else 

r™m _m 1 _ m-li 

X x k , — \x k , X fc _ 1 } 

end if 

Refinement 

Propose {at^} ~ q 2 {x k \x^,x\ 

Compute 4> 2 (u, x* k ,x™ 

=17T lo S 


n \ 
1/ 


u P{x^\x^_T)qi{xl\x^ ,x^_ x ) 

p(x* k \x^_ 1 )q 2 (x™\ xl ,x™_fi) 


Compute A 2 m ’ k (x™, x k ) and {pi(xf, x%)}^ with 
the routine described by Algorithm [3] 


if A 0 2 - k (x* kl xZ)>Mu,xlx™,a 

{*n = Ki 

end if 
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then 


if m > Nb V 
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(7TI Nburn ) 


then 

_ ryTT 
- "L h. 


end if 
end for 
end for 


p{x k \z 1:k ) = -A- K x k - 
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where Px(:) and pc (■) represent the likelihood of the target 
and clutter measurements respectively. Each individual mea¬ 
surement represents a point in the two dimensional observation 

r -1 T 


space, z k = 


^x,k’ * y,k 


. In the case of a measurement 


from the target, the likelihood is modelled as px{z k \x k ) = 
Af(z k ;x k: Y,). The clutter measurements are independent of 
the state of the target and are uniformly distributed in the 
visible region of the sensor, resulting in the clutter likelihood 
taking the form of pc(4) = u RS z i, k ) u R y { z i, k )- 




































Algorithm 3 Adaptive Subsampling Routine 

1: Given: The current and proposed states of the Markov 
chain, {x k , x k }, the complete measurement set, Zk = 


r„i ~ M k 
\ z k’ ' 


}, 5, and ?/>(•). 


2: Initialise: Number of sub-sampled measurements, 
S m ,k = 0, Approximate log likelihood ratio subtracted 
by proxy, A = 0, set of sub-sampled measurements, 

z k = 0, initial batchsize, 6 = 1, while loop counter, 
w = 0. 

3: Compute an upper bound for the range, R k , according 
to dTH] ). 

Compute the proxy, {pi(x k ,xl)}f^\, according to 
DONE = FALSE 
while DONE == FALSE do 


w — w + 1 
f S'm.fe + l,* 
l z k > 


b,* i 


".V}~ 

z* k =z* k U{zt’ k+1 '*,~ 
A=i(5 ro , fc A+E i=Sm , fc+ i 

Sm,k = b 
6 W = 06 


u/repl. z k \ z fc 
b ,*-1 


r p( 




Compute c according to (O utilising 5 W . 

b = "ySm,k M k 

if l A + TR Si =1 Pi(nk,x* k ) - > c or S m>k == 

Mk then 
DONE = TRUE 
end if 
end while 

return A and {pi(x k ,x* k )}™\ 


The Taylor approximations used by the proxy in ( IT5l > are 
dependent on the gradient and Hessian of the log likelihood 
for individual measurements. Substituting the terms for the 
target and clutter likelihood in (l20l> and taking the logarithm 
results in the log likelihood for each measurement having the 
form 

^{xk) = log ^A a -A S(z' k ; x k , £) + , (21) 

where Ac = R x x R y represents the clutter area. The gradient 
can then be computed as 

w = Aa-ST 1 ^ -x k )U(zl,x k ,^) 

\ X jV(zl,X k ,V) + %; ’ 

and the Hessian is given by 

-Aa£-W(z* ; x k , S) (a a -A f(zl; x k , S) + 0) 

H = - —2- h 

(AaA+ 0) 

- Xk ) (V-\zl - x k )) T Af(zU x k , S) 

(Aa -M{zl-x k ,T,) + As.) 

(23) 


B. Implementation Considerations 

The primary difference between the standard and adaptive 
subsampling sequential MCMC is that the latter requires less 
evaluations of the log likelihood. However, there are also 
additional computations which are introduced to achieve this. 
These calculations are minimal and typically performed for a 
fraction of the time spent on the calculation of the likelihood, 
when M k is sufficiently large, and are thus considered negli¬ 
gible. In this section we discuss these computations in more 
detail. 

The proxy, given in ( fl5] > is extremely efficient to compute 
in comparison to the log likelihood. This is conditioned on 
the availability of the gradient of the log likelihood (i.e. (l22l >I 
evaluated at a specific point. Currently, we only update this 
twice per time step (represented by line [5] in Algorithm [2]). 
Once, at the beginning of a time step, where the specific point 
used is the predicted mean of the Markov chain at the previous 
time step. Secondly, the current state of the Markov chain after 
the burn in period. As the number of MCMC particles, N, is 
typically several magnitudes larger than 2, these calculations 
are considered negligible. 

The calculation of an upper bound on the range in ( ITSl ) 
is also extremely efficient to compute in comparison to the 
log likelihood. This is conditioned on the availability of the 
maximum of the Hessian in (l23l) . In our application we found 
that the maximum of the Hessian is independent of the mea¬ 
surements and can hence be determined prior to the running 
of the algorithm (represented by line [2] in Algorithm [2]). 

The proposal distribution used for the joint draw step in the 
tracking scenario is defined as: 

l Np 

qi = p{x k \xk-i)—^2s(xk-i -x^i). 

P 3 = 1 

(24) 

The proposal distribution used for the refinement step in the 
tracking scenario is defined as: 

<?2 (®fc|®k l I *k l -i) =N(Xk',’Zq)- (25) 

The refinement step represents a local move. 


IV. Results 

Consider the scenario of a target moving through a highly 
cluttered environment. A sensor monitoring the target returns 
multiple target and clutter measurements at each time step. 
We applied the standard and adaptive subsampling sequential 
MCMC algorithms for the inference of the latent states of the 
target over several experiments with different parameters. 

Two different metrics are used to compare the performance 
of the algorithms. Firstly, the root mean square error (RMSE) 
of the position. The RMSE for each time step is calculated 
over a number of independent simulation runs according to 


RMSE = 



£(L 


2—1 


Xi) 2 , 


(26) 

















where X t represents the ground truth, X, represents the 
algorithm estimate, which corresponds to the mean of the 
N MCMC samples in this application, and Nj represents 
the number of independent runs. The RMSE of the states 
corresponding to the position are averaged to obtain a single 
result. The RMSE of the position illustrates the tracking 
accuracy of the two algorithms. 

The second metric is the normalized number of sub-sampled 
measurements required for likelihood calculations. 

JJ _ J_ (Sm,k)R (21} 

~ T 2 NM k 1 ’ 

k =1 

where {S m ^)jD and (Sm^R refer to the number of sub¬ 
sampled measurements from the joint draw step and refinement 
step respectively. The standard sequential MCMC algorithm 
requires to evaluate the likelihood 2 NM^ times at each time 
step, this corresponds to D = 1. Thus the D value is 
only shown for the adaptive subsampling sequential MCMC 
algorithm. It illustrates the fraction of likelihood evaluations 
which are required at each time step versus the standard 
sequential MCMC algorithm. 

A. Parameters 

The following parameters, unless otherwise specified, were 
used for all experiments. Simulation parameters: N = 500, 
N burn = 125, T = 20, Ni = 50, S, = 0.017. Motion 
model parameters: T s = 1, cr x = 0.5. Target observation 
model parameters: Ax = 500, £ = I. Clutter parameters: 
Ac = 2000, A c = 4 x 10 4 . Subsampling parameters: 7 = 1.2, 

6 = 0.1, p = 2. 

B. Performance Evaluation 

The first experiment illustrates the performance of the 
algorithms for different values of the mean total number 
of measurements in Fig. 1. The ratio between the mean 
number of clutter measurements and mean number of target 
measurements is fixed at 4:1. The RMSEs of the algorithms 
are in agreement, however, it is noted that an increase in the 
total mean number of measurements results in substantial com¬ 
putational savings. The amount of computational saving is as 
high as 80% with no significant loss in tracking performance. 

In Fig. 2 the ratio between the mean number of clutter mea¬ 
surements and mean number of target measurements is varied. 
This allows for the observation of the performance when there 
is a varied amount of information about the target present in the 
measurements. The RMSE results show agreement between 
the two algorithms with an increase in computational savings 
when the mean number of target measurements is higher. 

Fig. 3 illustrates the influence of varying the covariance 
matrix of the target observation model. The RMSEs of the 
two algorithms are in agreement. It is noted that a smaller 
computational saving is incurred as the measurement model 
becomes more precise. This result seems counter-intuitive. The 
reason for this is due to the Taylor approximation for the proxy. 
The upper bound for the range, , becomes a weaker bound 
as the observation model becomes more peaked. 



(a) RMSE comparison for different values of mean number of total 
measurements. The dotted lines represent the results from the standard 
sequential MCMC, and the full lines represent the results from the 
adaptive subsampling sequential MCMC. 



(b) Comparison of the normalized number of subsampled measure¬ 
ments evaluated in the adaptive subsampling sequential MCMC for 
different values of mean number of total measurements. 


Fig. 1: Performance comparison for a different mean 
number of total measurements with a constant clutter to 
target measurement ratio of 4:1. 

V. Conclusion 

In this paper, we presented an adaptive subsampling sequen¬ 
tial MCMC algorithm for target tracking. We have shown that 
this approach results in substantial computational savings when 
there is a large number of measurements, and most importantly, 
without sacrificing tracking performance. 

There is a wide scope for future work. From an application 
perspective, considering a multi-target scenario with different 
levels of clutter, and also from an algorithmic perspective, 
further research on the influence and implementation of a more 
efficient proxy. 





















(a) RMSE comparison of different values of mean target measure¬ 
ments. The dotted lines represent the results from the standard 
sequential MCMC, and the full lines represent the results from the 
adaptive subsampling sequential MCMC. 



(b) Comparison of the normalized number of subsampled measure¬ 
ments evaluated in the adaptive subsampling sequential MCMC for 
different values of mean target measurements. 


(a) RMSE comparison for different measurement covariance matrices. 
The dotted lines represent the results from the standard sequential 
MCMC, and the full lines represent the results from the adaptive 
subsampling sequential MCMC. 
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(b) Comparison of the normalized number of subsampled measure¬ 
ments evaluated in the adaptive subsampling sequential MCMC for 
different covariance matrices. 


Fig. 2: Performance comparison for a different number Fig. 3: Performance comparison for different covariance 

of mean clutter to target measurements ratios. matrices where S = cr?/. 
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